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ABSTRACT 

We discuss B c mass (IS state) and decay constant fs c calculated by lattice non- 
relativistic quantum chromodynamcs(NRQCD) method. In leading order of v 2 } we found 
that Mb c = 6.33(2) GeV and fs c = 395(2) MeV where the error bar is the statistical error 
only. Using these values, we estimate QCD effects to leptonic decay width of B c . The decay 
width is given by T(B C l + u x ) = 0.86(15)mf (1 - 0.0250(2)mf f x I0" 14 (GeV), where m, is 
in GeV. Relativistic correction and a s correction to fs c have also been considered. 



Theoretical understanding on hadronic and electromagnetic annihilation decay of 
heavy quarkonium system gained a solid footing due to a new factorization theorem by 
Bodwin, Braaten, and Lepage [1](BBL). They observed that in heavy quarkonium decay 
processes, typical short distance scale, M~ l (Compton wavelength of heavy quark), is well 
separated from long distance scale, (Mv)~ l (size of heavy meson) due to small quark ve- 
locity v. Then using velocity scaling law of various operators in non-relativistic QCD[2], 
they showed that decay rates can be written as a sum of factorized products of perturbative 
coefficients and non-perturbative matrix elements. Compared to earlier calculations, their 
argument is more consistent because they showed how to do systematic expansion in v 2 and 
a s which makes the factorization work in every order. It is valid even for P-wave case (or for 
higher orbital states), where explicit perturbative calculations questioned naive expectation 
by showing an infra-red divergence in (9(a^)[3]. 

In this work, we are interested in various matrix elements involved in BBL factoriza- 
tion theorem. The matrix elements in this factorization contain non-perturbative information. 
Such information may be obtained from either fitting experimental data or doing lattice QCD 
simulation. When there is no available experimental data for the matrix elements, predictions 
by lattice QCD simulation will be useful. Weak annihilation decay of charmed B meson, B C} 
is one of such heavy quarkonium systems to which the new factorization can be applied and 
on which has no experimental data yet. Although there exist studies of B c system based 

on models [4, 5, 6], we think that the result from lattice calculation is worthwhile because 
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lattice calculation is based on first principles of field theory so that the result can be im- 
proved systematically in terms of accuracy if desired. Interestingly, there are expectations 
that B c is within experimental reach of LEP or Tevatron and most certainly of LHC. It is 
anticipated that decays of B c will reveal rich information on the Standard model because 
total annihilation decay of B c can proceed only through the weak interaction [7, 8, 9]. 

We have been calculating the matrix elements associated with heavy quarkonium decay 
such as cc and bb system using lattice formulation [10]. Extension of our calculational scheme 
to charmed B meson system is straightforward. Let us describe our method briefly in the 
following (further details on our lattice simulation method can be found elsewhere [12]). We 
use non-relativistic lattice QCD (NRQCD) formulation in which the scale separation for the 
new factorization theorem is most transparent. NRQCD Lagrangian upto v 2 } 



is employed in our calculation. Thus our result will have 0(v 4 ) corrections, which may be 
as large as 20 %. Also, we cannot address fine structure of the spectrum since such effects 



the quenched approximation which neglects light quark vacuum polarization effects. We use 
200 gauge field configurations on 8 3 X 32 lattice volume at (3 = 5.7 which were generated by 
use of Metropolis and over-relaxation algorithm. Although the spatial volume of our lattice 
configuration is small (~ 1.6 fm), we think that the finite volume effect in our calculation for 






means just gluon degrees of freedom since we work in 
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B c meson system will be small because the spatial size of B c itself will be small. We fixed 
the gauge field to Coulomb gauge prior to the propagator and the matrix element calculation. 
Under these quenched gauge field background, the quark propagator in lattice NRQCD is 
calculated as the following [2], 

G(x, t + l) q = (l- H /2n) n U\x } t)(l - H /2n) n G(x } t) q + ( 2 ) 

with G(x,t) q = for t < 0, 
where 

H = -A 2 /2M - 3(1 - u )/M . (3) 

A 2 is the lattice covariant discrete laplacian. The second term (= E su b) is the energy shift 
that arises from the deviation of the plaquette from 1 in mean field theory in "tadpole" 
improvement scheme [11] and u is defined by 

u = (0\^TrU plaq \0)^ (4) 

where U p i aq is the plaquette value. We choose n = 2, E su b = 0.605524909 for c quark, and 
E su b = 0.154739082 for b quark and u = 0.8607297. Bare c quark mass parameter is chosen 
to be 0.69, and that of b quark mass is 2.7 in lattice unit. These quark masses are used 
because they reproduce J/^ mass and T mass at (3 = 5.7. Actually, we need to adjust two 
quark masses so that the correct B c meson mass is reproduced. However, since there is no 
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experimental data yet, we use these quark masses. Then, B c propagator is given by, 

G(t) Bc = J2G(x,t) q G\x,t) q , (5) 

X 

Fig 1 is the behavior of B c propagator. We fit G(t)s c to the following form, 

G(t) Bc = Ae~ E \ (6) 

CERN Library MINUIT is used for minimization of the correlated \ 2 ■ We found that the 
error bar from the single elimination jackknife method is similar to that from the \ 2 method in 
magnitude. Fig 2 shows effective mass plot. Each points is logarithm of the ratio G(t) /G(t-\-l) 
and the error bar is from the jacknife method. The plateau value is 0.955(1) in lattice unit. 
B c mass is given by, 



M b - C , c i = \(M bl - E hl a~l + M c - C - E c - c a~J) + E^a^ (7) 



where 



M c - C = -Mj h ^ + ±M„ C . (9) 



r/b has not been seen but there exists a good estimate[13]. From PDG [14] and NRQCD group 
estimate, M cE = 3.06739(56) GeV and M hl = 9.45299(142) GeV. From our calculation[12], 
E h i = 0.7948(5), E C c = 1.014(2) in lattice unit. Here, a~^{ = 1.37GeV) is the lattice spacing 



for bottomonium, and a c ^ (= 1.23GeV) is that for charmonium. We estimate the lattice 
spacing for B c , a^ 1 or by 1.30 = ±(1.37 + 1.23) or 1.30 = ^1.37 x 1.23 in GeV. This 
ambiguity in scale is just due to the quenched approximation which can be remedied in the 
future simulation. The introduced error by this estimate of the lattice spacing is more or less 
the same in magnitude as the errors resulted from 0(v 4 ) terms which is not included in our 
calculation. Thus, by adopting a^"- 1 or a"* = 1-30 GeV as the scale in our results and by 
using the plateau value E h ^ c i = 0.955(1), we get M Bc = 6.33(2)(GeV). 

The decay constant, f Bc} is defined by 

f B y = (0\bY7 5 c\B c (p)). (10) 

In NRQCD, 

M B J Bc ~ C o (0\x\4>c\B c ) + C 2 (0\(D Xb y ■ D^ C \B C ), (11) 
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ith the normalization, (B c (p) \B c (p')) = (2ir) 3 2p S 3 (p—p'). The second term is the relativistic 



correction. The coefficients Co was calculated upto 0(a s ) in [9]. It is 



c ° = 1+ fi^m log 5s- 2 i- (12 > 



Also 



_ (M b + M c f 

62 " " 8(M bMc y ■ (13) 



With M 6 physical = 4.5(GeV),M c P h y sical = 1.5(GeV) and a s (fg£) = 0.34, C = 1.0(= 0.85) 



without a s correction (with a s correction) and C\ = — 0.09668(GeV 2 ). 

Therefore, the matrix elements which is needed for B c decay rates are, 

G 1 = (0\xtil>c\B c C S)) = yfiMZGi. (14) 

F x = (0\{D Xb y ■ D^ c \B c ( l S)) = yfiM^K (15) 
The non-covariant form of F\ is, 

K = (OKVx^) 1 • V*p c \B c CS)) = ^2M^F{- (16) 

After removing the exponential fall-off from the binding energy, the asymptotic value of 
B c propagator gives the wave function at the origin, ^^(O)) 2 = 0.1125(6)a~ 3 . From the 

2 — 2 

definition G\ = ^y|i?(0)| 2 , this translates into = 0.1125(6)a~ 3 . Fig 3 shows the ratio 

of ^p- to the B c propagator, where the plateau value gives 2F[/G\ = 1.813(2)a . On the 
other hand, we get = 2.400(2)a" 2 . 

In leading order of a s (or v 2 ) } consistency tells us that fs c = 395(2) (GeV) from 
f Bc = Beyond leading order, to be consistent in perturbation expansion, a s correction 

in Co should be considered, and the second matrix element (relativistic correction term) 
need to be included and the lagrangian should contain 0(v 4 ) terms because a s ~ v 2 [l]. 
In addition, MS matrix elements need to be extracted from the calculated lattice matrix 

elements since the perturbative coefficients, Co and Ci, in decay constant is calculated in 
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MS scheme [9]. This requires operator matching between these two different regularization 
scheme. Such operator matching is similar to bb system[15] except the fact that the mass of 

quark is different from that of anti-quark. For G^ s , the feynman diagram, Fig 4a, needed 
to be evaluated with zero external momentum and the scale of the diagram needed to be set. 
The algebraic expression for the diagram in coulomb gauge is given by 

= l*r f d3q 1 ( 1 - ^ (W\ 

where fi = l — ^-sin 2 y(M; is either charm quark (0.8016) or bottom quark (3.137) mass 
in lattice unit). The numerical evaluation of the integral gives — 0.1 28298 (4) a s . By use 
of a s = 0.188(p = ir / a) at (3 = 5.7, the integral is equal to —0.0241220(7) numerically. 
Since the modification due to matching is small, we do not do elaborate scale fixing for the 
diagrams. If we take the modification due to matching and a s corrected coefficient Co into 

consideration, we get fg^ = j^jj^- = 344(2) MeV from the first term of (11). Since our G\ is 

0(v 2 ) lagrangian result, the above modification is not directly useful in the phenomenological 
discussion. It just gives us a rough idea on the magnitude of each contributions. On the 
other hand, similar consideration on our F\ is relevant. 

For F^ s , we need to subtract a divergent piece which is proportional to v^?/>(a) 
from feynman diagrams, Fig 4b, 4c, and 4d. Such piece can be deduced from the external 

momentum independent portion of these diagrams. The continuum D-D 1 is the p-p' dependent 
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portion of these diagrams, which can be obtained by taking zero external momentum limit 
after taking derivative with respect to external momentum. Similarly, for F^ s \ a divergent 
piece from feynman diagrams, Fig 4b and 4c, need to be subtracted and the continuum V • V 
is the p ■ p' dependent portion of these diagrams. Thus 



F r = _2_ (f i__L_^ (18) 



where e is p ■ p' dependent portion of the diagram 4b and 4c, and / is the subtraction 
piece. With a s (ir/a) = 0.188, numerically, e = (2.6866(1) - 0.35689(7))a s = 0.43798(9) 

and / = (2.6743(1) + 5.1913(4))a s = 1.4787(1) for Ff and = 2.6743(l)a s = 0.50277(2) for 

Ff'. Using Ff\ we get = 0.324(2)(GeV) 7/2 . On the other hand, using F~£, we get 

— 0.260(10)(GeV) 7/ ' 2 . If we use the non-covariant operator result, relativistic correction gives 

(2) (2) p 

fg> = -17.6(l)MeV from fg> = C 2 ^. However, since the difference between the covariant 

result and the non-covariant result is significant (even the sign is different), extracting F^ s 
from lattice study needs further study. Thus in the following discussion, we use only leading 
order result. 

By use of the leading order Mb c and /b c , the leptonic decay width [9], 

1 Til 2 

Y(B C -+ l+ut) = -\V hc \ 2 GlM B Jlm*(l - ) 2 , (19) 



becomes 

T(B C -+ l + vi) = 0.86(15)mf (1 - 0.0250(2)m?) 2 x 10" 14 (GeV), (20) 

where V bc = 0.040(5), Gp = 1.16639(2) x 10- 5 GeV -2 [14] and m, in GeV have been used. 
Since the major source of errors in the calculated decay width is V& c , we will get better 
information on V bc once the leptonic decay width is measured experimentally. 

In conclusion, using non-relativistic lattice formulation of QCD, we calculated IS 
state B c mass (= 6.33(2)GeV) and the decay constant /b c (= 395(2)MeV). Using these 
non-perturbative information, the leptonic decay width of B c has been discussed. 
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quantity 


our result 


Eichten and Quigg 


Chang and Chen 


Kiselev et al 


M Bc (GeV) 
f Bc (MeV) 


6.33(2) 
395(2) 


6.194 ~ 6.292 
479 ~ 687 


6.34 ~ 6.35 

422 ~ 450 


6.301 ~ 6.344 
456 ~ 510 



Table 1: comparison of Mb c and fs c from various studies 
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Figure captions 



1. B c propagator as a function of t in log scale. 



2. effective E h7 = function of t. 



3. ratio of the matrix element 2F[ to B c propagator as a function of t. 



4. feynman diagrams necessary for operator matching. Vertex factor for a) and d) are 

unit operator, and that for b) and c) are lattice V • V. Dotted line means temporal 
component of the gauge field and wavy line means spatial component of the gauge field. 
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